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Our concern is the behaviour of the elementary cellular automata with state set {0, 1} over the cell set 
Z/nZ (one-dimensional finite wrap-around case), under all possible temporal rules (asynchronicity). 

Over the torus Z/nZ (n < 10), we will see that the ECA with Wolfram update rule 57 maps any 
v e Fj to any w S FJJ, varying the temporal rule. 

We furthermore show that all even (element of the alternating group) bijective functions on the 
set FJj = {0, . . . , 2" — 1 }, can be computed by ECA-57, by iterating it a sufficient number of times with 
varying temporal rules, at least for n < 10. We characterize the non-bijective functions computable 
by asynchronous rules. 

The thread of all this is a novel paradigm: 

The algorithm is neither hard-wired (in the ECA), nor in the program or data (initial configura- 
tion), but in the temporal order of updating cells, and temporal order is pattern-universal. 
Keywords: Cellular automata, asynchronous, update rule, universality. 

1 Introduction and Notation, Asynchronicity 

We consider elementary cellular automata, i.e. with state set 5 = F2 = {0, 1} and update neighborhood 

(d-i, Ci,Ci + i) for cell c,. 

The cell index (site) i will come from Z/raZ for some n > 3, i.e. we consider the finite one- 
dimensional torus, indices wrap around. In Section 2, we consider patterns "How universal can a map- 
ping on ¥2 become?", and Section 3 covers functions FJJ B v — > w G F!>. 

The 256 ECA's group into 88 classes under the symmetries 0/1 and left/right neighbor, see Ap- 
pendix A. It is sufficient to consider one member per class. 

The Wolfram rule ECA = Yl=o 2 ' Pk e {0> • • • > 255} defines the behavior. A cell with neighborhood 
(c;_i,c,,c,_|_i) G F|, summing up to k := 4c ; -_i +2c,- + c, + i £ {0, ... ,7} is replaced by := pk- 

Example 1 The behaviour of the ECA with Wolfram rule 57io = 00111001 2 is given in Table 1. We 
have that 0c,T h-> Cj, all other cases 0c ( 0, lc,0, lc,T H> c;. 



Ill ^0, 
110^0, 
101 ^1, 
lOOi-^l, 



011 h-H, 

010 ^0, 

001 ^ 0, 
000 1 — y 1. 



Table 1: ECA-57 
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1.1 State-of-the-Art 

The study of asynchronous cellular automata started with Ingerson and Buvel's 1984 paper Q. 

Lee et al. give an asynchronous CA on the two-dimensional grid ZxZ, which is Turing-universal. 

Fates et al. [1] consider ECA's with quiescent states (000 1 — )• 0,111 i — >• 1 , i.e. with even Wolfram rule 
> 128). They consider fully randomized ECA's. 

A good overview is given in the thesis H by Sharkar. 

Nevertheless, all these articles treat asynchronous CAs with randomized clocks. 

Our concern is instead the (fully deterministic) behavior of a suitable ECA, with any fixed initial 
configuration, under all possible temporal sequences. There seems to be no work on the combined effect 
of all deterministic temporal rules, synchronous and asynchronous, so far. 

Definition 1. Temporal Rules — Asynchronicity Rules 

Let the set AS,, of asynchronicity rules over Z/nZ consist of all words of length n over the alphabet 
{<,=, >} such that both < and > occur at least once. We also include the word "= • • =", the syn- 
chronous case, and have AS,, = ({<,=, >}"\({<,=}" U {=,>}")) U {="} with \AS n \ = 3"-2 M+1 +2. 

Given a rule AS = ASo • • • AS„_i, AS, = "<", "=", and ">", resp., defines that cell c; updates after, 
simultaneously with, resp. before c i+ i, V0 < i < n — 2. AS„_i refers to cell c„_i with respect to cq. 

For any partition (Si,.. -S m ) of the cell sites, i.e. 0^ =l Sk = {0, ...,n — 1}, let its temporal rule be 
AS; = <,=, > , resp., if i G S l rf j , i+ 1 G Si(;+i), and i(i) is >,=, <, resp., than i(i + l) (we say that site 
i is "bigger" if it comes before i + 1, hence dominates it). 

With the exception of =" (synchronous case), both < and > must occur at least once, since otherwise, 
by wrapping-around, each cell would update only after itself and the temporal rule would thus not be 
well-defined, e.g. =<= leads to c\ with C2 after C3 with c\, so c\ after, and thus before, itself. 

Example 2 Let n = 4, and AS = "<=>>": Cell updates after cell 1, 1 with 2, 2 before 3, and 
3 before 0. Hence the temporal order is (1,2|3|0), first 1 and 2 simultaneously, then 3, finally cell 0, 
i.e. Si = {1,2}, 52 = {3}, S3 = {0}. Analogously, "<><>" leads to (1,3|0,2), and ">=><" leads to 
(0|1,2|3). 

One might be inclined to partition the n cells into sets Si,... ,S m C Z/nZ, and update those in Si first, 
then cells from S2 and so forth. This, however, is too fine-grained: 
Theorem 1 

Consider two partitions (Si,. . . ,S m ) and (S[, . . . ,S' ,) of the cell set Z/nZ and define functions 
a,b,c,a' ,b' ,c' such that 

V/ G {0, ... ,7i — 1} eS a (i),ieS b ^,i + l eS c ^,... 

■ ■ ■ i ~ 1 e s 'a'(i) > i e s 'b'({) > i + 1 G s 'c'(i) ■ 

Then, ifsga(a(i) —b(i)) = sgn(a'(/) — b'(i)) and sgn(/3(i) — c(i)) = sga(b'(i) — c'(j)),Vi, i.e. the relative 
update order of cells i — 1 , i, i + 1 is the same for S and S', then updating according to S or according to 
S' leads to the same result, and this is described by the following asynchronicity rule (Table 2). 

Proof. By construction. Since the relative temporal order of cell c; with respect to c,_i and c,+i is the 
same for (St) and (S' k ) by sgn(a(/) — b(i)) = sgn(a'(/) — b'(i)) and sgn(b(i) —c(i)) = sga(b'(i) —c'(i)), 
both partitions lead to the same overall behaviour, which is described by AS. □ 

The construction by the theorem shows that the AS G AS,, are sufficient to distinguish the behaviour. 
On the other hand, all these AS are necessary and can lead to different behaviour (at least for some 
ECA's), since any AS; 7^ AS'; will lead to a different order of updating cells c; and c ;+ i. 
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sgn(a-ft) 


sgn(b — c) 


AS,-_i 


AS; 




-1 


-1 


> 


> 


a before b before c 


-1 





> 




a before b with c 


-1 


+ 1 


> 


- 


a and c before b 





-1 


= 


> 


a with b before c 












a with b with c 





+1 




< 


c before a with b 


+1 


-1 


< 


> 


a before b and c 


+1 





< 




b with c before a 


+1 


+1 


< 


< 


c before b before a 



Table 2: Local asynchronicity 



Example 3 For n = 6, "<><><>" requires the odd cells 1,3,5 to update before the even ones 
0,2,4. There are 13 partitions of three elements, e.g. (1,3,5), (1|3, 5), (1,5|3), and (5|1|3), and thus 
13 2 = 169 partitions (S k ) for this AS. 

Definition 2. By ECAas (v) = w, we mean that the elementary CA with rule EC A maps v G {0, 1}" 
to w G {0, 1}" via the temporal sequence AS. 

Example 4 ECA-57 <><> (1000) = 1110, in two steps: 1000 ^ 1100 ^ 1110, where underlined 
cells are active in the next step. 



2 The Finite Torus Z /nZ: Patterns 

In this section, we work on the torus Z/nZ, and consider all EC As for all initial configurations. We 
apply a fixed temporal rule AS G AS n repeatedly, t times, and ask, whether these 5 pattern universality 
properties hold: 



(o) 

(0 

(ii) 

(Hi) 
(iv) 



3vGF^, VwGF!J, 

VvGF^, VwGF^, 

Vv G FJ, 3t G N, 

3t G N, Vv G F%, 

3t G N, Vt > To, 



3tgN,... 

3tgN,... 

VwGF^,... 

VwGF^,... 

Vv,wGF",.. 



> 3 AS G AS n ■ ECA^ S (v) = w. 



All results are experimental i.e. derived from exhaustive computer simulations for the stated lengths. 
We start with 

(o) 3v G F£, Vw G F^, 3t G N, 3 AS G AS n : ECA^ S (v) = w. That is from some v we eventually reach 
any w. We give the largest number of w's reached for some v, for n = 4, 8, and 12. To satisfy (o), these 
must be (16,256,4096). 

The 3 EC A families (1,1,1), 200 (1,1,1), and 204 (1,1,1) are resilient to asynchronicity. They have a 
constant result, for all n. 



ECA-51 (2,2,2) varies between at most two results. 
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The next 49 EC A families are ordered by increasing image size for n = 12: 



140 


(2,6,16), 


160 


(12,130,1182), 


164 


(13,197,2930), 


108 


(16,256,4052), 


136 


(2,9,27), 


2 


(11,211,1477), 


24 


(15,211,2961), 


56 


(16,256,4066), 


128 


(2,16,49), 


72 


(11,131,1499), 


34 


(13,209,2998), 


74 


(15,255,4071), 


132 


(4,18,81), 


76 


(11,131,1499), 


130 


(14,211,3160), 


73 


(16,256,4084), 


32 


(7,31,127), 


172 


(11,137,1506), 


94 


(16,216,3448), 


33 


(16,256,4092), 


8 


(5,45,320), 


168 


(12,147,1601), 


152 


(14,237,3561), 


10 


(13,253,4093), 


4 


(7,47,322), 


13 


(16,168,1792), 


138 


(13,238,3751), 


134 


(15,255,4093), 


12 


(7,47,322), 


232 


(12,156,1830), 


104 


(14,232,3824), 


42 


(15,255,4093), 


28 


(11,91,641), 


77 


(12,156,1830), 


162 


(16,250,3970), 


35 


(16,256,4094), 


29 


(12,92,642), 


142 


(12,140,1848), 


170 


(16,256,3976), 


43 


(16,256,4094), 


44 


(12,100,870), 


78 


(15,167,1851), 


15 


(16,256,3976), 






156 


(4,64,1024), 


36 


(14,162,1943), 


150 


(12,240,4032), 






40 


(11,119,1052), 


5 


(16,216,2542), 


1 


(16,256,4051), 







The 4 EC A families 6, 14, 18 [for n > 7], 50 [for n > 4], miss exactly one pattern, leading to 2" — 1 
in general. 

Finally, the 31 EC A families 
3, 7, 9, 11, 19, 22, 23, 25, 26, 27, 30, 37, 38, 41, 45, 46, 54, 57, 58, 60, 62, 90, 
105 [?2 /0 mod 4], 106, 110, 122, 126, 146, 154, 178 [n^3], 184 [n^3], 
satisfy property (o) (for 3 < n < 12). 

(?) Vv,w G F^ElAS G AS n ,3x G N : ECA^ s (v) = w. From the 31 families satisfying (o), most fall 
short for some v. We give the smallest number of w reachable from some v, for n = 4,8, and 12, this 
should be (16,256,4096) to satisfy (?). 

Eighteen ECA families are insensitive (or resilient) to asynchronicity for at least some v, the same w 
resulting for all AS. Hence, (1,1,1) patterns are reached: 
22, 26, 30, 38, 46, 54, 58, 60, 62, 90, 106, 110, 122, 126, 146, 154, 178, 184. 
ECA family 7 reaches 2" — 1 for n ^ mod 3 and only 1 pattern for n = mod 3. 
ECA family 45 has 2" — 1 patterns for odd n, 1 for even n. 

Six ECA families get near the full 2" for all w: 3 (15,233,3411), 9(12,243,3963), 11 (15,233,3515), 25 

(16,251,4031), 27 (16,253,4052), 43 (12,236,3554). 

The following 6 ECA families satisfy (?) at least for certain [n] (3 < n < 12 considered): 

19 [3-12], 23 [3,5,7,9,11], 37 [4-5,7-8,10-11], 41[3,5,7-12], 57[3-12], 105 [3,5-7,9-11] all generate 

2" patterns for these [n] . 

(ii) — (iv) From now on, we will consider the 6 ECA families satisfying (?'): 19, 23 (n ^ mod 2), 
37(«^0 mod 3), 41, 57, and 105 («^0 mod 4). 

(?'?') Vv G ¥1,3t G N,Vw G F^Bas G AS n : ECA^ s (v) = w; i.e. for fixed v, all w are reached at the 
same time. 

We considered % up to 20000, and obtain: 
ECA-23: No v has any T < 20000 to satisfy (??). 
ECA-19,-37,-41: For some v, there is no % < 20000 to satisfy (??). 

ECA-57 satisfies (??), for n > 5 and all v. The largest T required is 28 for n = 5; 14 for n = 6; 10 for 
7 < n < 13; and 9 for n = 14 and 15. 

ECA-105 satisfies (?'?') for odd n > 1 and all v. The largest x required is 30 for n = 7; 16 for n = 9, 1 1 , 13; 
and 8 for n = 15. 

In general, the time x decreases with n, since the number of patterns, 2 n , increases slower than the 
number of asynchronicities, 3"-2 n+1 +l, and thus for larger n, AS n is more likely to satisfy (?'?') early on. 
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(Hi) 3t G N,Vv,w G F^Bas G AS„ : ECA AS (v) = w. All transductions v i->- w are done in the same 
time. 

From the result of (ii), we can infer that at most ECA-57 and ECA-105 can satisfy (Hi). 
ECA-57 has a joint T„ at which all transductions vh-w are satisfied simultaneously in these cases: 
T 5 = 445, T 7 = 70, T 8 = 242, t 9 = 35, no = ■ ■ ■ = t m = 13, Ti 5 = 10. 
For ECA-105, we have T 7 = 570, T 9 = 14, Tn = T13 = 6, and T15 = 8. 

(iv) 3t G N,Vt > T ,Vv,w G F^,3as G AS n : ECA^ s (v) = w. Eventually all transductions vh-w 
can be done at all times. 

Theorem 2 

There is no To G N, Vt > To, Vv,w G F!>, 3 AS G AS n : ECA AS (v) = w, i.e. (iv) can not be satisfied. 
Proof. Consider the case w = v. 

For each rule AS, applying AS repeatedly, starting at v, we will either return to v at some time, which 
is the period length per(AS), the length of the cycle of AS containing v, or else v is in a preperiod and 
will never be reached again. Therefore, either v is in the preperiod and thus will not reappear, or else 
there is no preperiod, and v = w appears exactly after k ■ per(AS),\/k, applications of AS. 

Let now PER = lcm^$(per(KS)), where AS runs over those temporal rules without preperiod. Ap- 
parently, after k ■ PER(v),VA:, applications of AS, we return to v, for all these rules without preperiod 
simultaneously. After k ■ PER(v) ± 1 steps, \/k, we are not at v (unless the period is 1, and thus v is a fixed 
point). Hence, v i-> v is impossible for all these timesteps k ■ PER(v) ± 1, and there is no such To- 

Finally, in the case that v is a fixed point under AS, no cell changes its contents for this temporal rule 
and thus only w = v, but no w 7^ v is ever reached. □ 

Definition 3. We call an elementary cellular automaton pattern-universal at length n, if it is able 
to convert any pattern v in into any other, i.e. satisfies property (?) (Vv,w G FI}, 3as G AS n , 3t G N : 
ECAX s (v) = w). 

If an ECA is pattern-universal for all n > no, for some no, it is called uniformly pattern-universal. 
Corollary ECAs from the 6 families 19,23,37,41,57, and 105 are pattern-universal for the lengths 
indicated in property (i) above. 

We conjecture that ECAs from families 19, 41, and 57 are uniformly pattern-universal. 

3 The Finite Torus Z /nZ: Functions 

In Section 2, we focussed on transductions vHw, which — in general — used different temporal rules 
for different v's and w's, but for each pair (v, w) stayed with the same rule, applied repeatedly. 

In this section, we are interested in functions F2 3 v >->■ f(v) = w G FJj, which use the same temporal 
rule sequence for all v, but — necessary to generate enough variation — concatenate several different 
temporal rules. 

We may restrict ourselves to ECA families 19, 23, 37, 41, 57, and 105. Given a function on F!J 
defined by the values w(v) G F!>, Vv, our question is thus: 

3k G N, 3ASi, . . . AS* € AS n , Vv G F^ : ECA AS ^(. . . (ECA AS[ (v)) . . . ) = w(v)? 
3.1 Bijective Functions 

We first consider bijective functions on FJ(. In this case the equivalent group-theoretic statement is: 

Do the ECAas G AS n (written as permutations on the set {0, 1, ... ,2" — 1}) generate the full sym- 
metric group S2" ? 
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To answer this question, we used the program GAP (Graphs, Algorithms, Programming) from RWTH 
Aachen (Prof. Neubiiser's group) and St. Andrews University 0. Thank you! 

We ran GAP on some subsets of only 3 asynchronicity rules to show that ECAas £ AS n generates at 
least the alternating group A2*, for 3 < n < 11. 

Trying directly to obtain the group generated by the full set ECAas £ AS n overburdens GAP already 
from n = 4 on. Therefore, in order to check for the generation of S2", it is then sufficient to exhibit at 
least one odd permutation, which is the case for n = 3, with the whole 5 2 3 generated — or to show that 
all permutations generated by ECAas £ AS n are even, which is the case for 4 < n < 11, and thus only 
Ay, but not S2", is generated in these cases. 

Out of the 6 ECA families satisfying property (/), ECA-57 and ECA-105 are the only ones, which 
have a locally bijective update rule. Therefore, only these families must be considered. We immediately 
have that temporal rules avoiding the symbol "=" are bijective, when the temporal rule is bijective, since 
different applications of that temporal rule do not interfere with each other. On the other hand, for n ^ 3, 
all temporal rules involving the symbol "=" lead to non-bijective functions, see next subsection. 

The rules excluding = define bijective functions, whenever the ECA itself is (locally) bijective, that is 
the application of such an temporal rule for a single cell yields bijectivity. Those temporal rules including 
= define the non-bijective functions. Hence, the only way to generate bijective functions for n > 4 is by 
using ECA-57 or ECA-105, and only applying temporal rules from {<, >}". 

ECA-57: GAP tells us that the 2" - 2 temporal rules from {<, >}"\{<", >"} always yield at least 
the alternating group Ay, which is in fact generated already by 3 of the temporal rules, for 3 < n < 10. 

The case S2" vs. Ay is easiest checked by hand: Is there some odd permutation within the temporal 
rules? This is only the case for n = 3. For 4 < n < 10, all temporal rules yield even permutations and 
thus can not generate the full Sy ■ 

Hence, for n = 3, all bijective functions are generated through ECA-57 by concatenation of suitable 
temporal rules, while for n > 4, only the even permutations from Ay (that is half of the 2"! bijective 
functions) are generated. 

ECA-105: GAP tells us that all bijective temporal rules combined generate only fairly small groups: 
< AS 3 > has order 24, | < AS 4 > | = 48 = 4! • 2 1 , | < AS 5 > | = 1920 = 5! • 2 4 , | < AS 6 > | = 11520 = 
6! • 2 6 , I < ASj > I = 322560 = 7! • 2 6 , all are far below |52»| =2"!, the number of bijective functions. 

3.2 Non-Bijective Functions 

We now turn to non-bijective functions. Then Im(f) C F'-J with \Im(f)\ strictly less than 2". 

We start with n = 3. The convex hull over all AS G AS3 has cardinality at least 2 3 !/2 = 20160 for 
the following ECAs, Table 3 (the other ECA with bijective update rule, ECA-105, generates only 344 
functions): 



ECA-25: 

ECA- 110: 

ECA-30: 

ECA-3: 

ECA-57: 

ECA-11: 

ECA-62: 



22496 
23166 
25258 
39155 
40320 
52934 
62683 



ECA-46 
ECA-41 
ECA-38 
ECA-27 
ECA-35 
ECA-54 
ECA- 19 



89110 
210493 
223102 
268034 
751760 
1.190.449 
3.519.992 



Table 3 



Image size for ECAs on 
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There are 8 8 , about 16 Mio., functions on F|. Hence, for n = 3, none of the ECA's even generates 
a quarter of all functions. The case ECA-57 is special in that this ECA actually generates all bijective 
functions, but no non-bijective one, for n = 3. 

In the sequel, n > 4, we consider only ECA-57, which has sufficiently many bijective functions, 
namely 2 n !/2, at least for 4 < n < 10. We will generate a considerable subset of all functions by suitably 
interleaving bijective and non-bijective temporal rules for ECA-57. 

We now consider ECA-57 for a temporal rule with a single = on Z/raZ, n > 4. 

Considering larger neighborhoods, with 2 cells changing simultaneously, also ECA-57 becomes non- 
surjective (we show the effect of ASi = "=" on the two middle cells for all configurations of 4 adjacent 
cells): 



v ^ ECA-57 (v) 



0000 




0110 


0001 




0101 


1000 




1110 


1001 


!->■ 


1101 


0010 


!->■ 


0000 


0011 


H> 


0011 


1010 


1 y 


1100 


1011 


I-)- 


1111 


0100 


!->■ 


0010 


0101 


!->■ 


0011 


1100 


1 y 


1010 


1101 




1011 


0110 


!->■ 


0100 


0111 


!->■ 


0101 


1110 


!-)■ 


1000 


1111 


1 y 


1001 



Table 4: ECA-57: Effect of ASi = "=" 

We obtain the patterns 0011 and 0101 twice, while missing 0001 and 0111. Hence the image is 
smaller than the full 2 4 by 2, or by a factor of 7/8. 

Extending this neighborhood of = to any size n, and using only < and > for the other n — 1 positions, 
before and after the = transition, ECA-57 behaves bijectively. Therefore, the whole image shrinks by 
just the factor 7/8, when applying = once. 

Since all temporal rules without = are bijective, and inclusion of more than one = shrinks the image 
even further, we have the following result on the functions that can be represented by ECA-57: 

Theorem 3 

Let the patterns from {<, >}"\{<", >"} generate at least the alternating group (which is the 
case at least for 3 < n < 10). 

Let f : ¥2 — > FJJ, n > 4 be any non-bijective function on at least 4 symbols. Let #(w) = | {v|/(v) = w} \ 
be the number of configurations v leading to configuration w. Then f is representable by ECA-57 under 
asynchronicity, if and only if 

£ L#( W )/2j>2- 3 . 
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Proof. We first introduce the functions # on F!J and @ on No: 

The multiplicity #(w) tells us, how often w is reached, i.e. is the size of the preimage of {w}. 

For k G No, let G No be the number of results w appearing with multiplicity k, @(k) = \ {w G 
¥2 : = k}\. In particular, @(0) = 2" — \Im(f)\ is the number of words avoided by the image of /. 
We have ^ k k-@(k) = 2 n . 

We make use of the temporal rule AS* := "<=> ••• >" which maps 2"~ 3 pairs (vi,V2) onto 2"~ 3 
words w, and otherwise is 1-to-l. Hence, for AS*, we have @(1) = 6 • 2 n ~ 3 , @(0) = @(2) = 2"~ 3 . We 
generate / by a chain / = 7Tj : oAS*o---o7r2 AS*o7z: 1 oAS*, alternating AS* and permutations n k G 52". 

For the second and every further application of AS*, we will join 2"~ 3 — 1 words vi with #(vi) > to 
2"~ 3 — 1 words V2 with #(v2) = 0, hence without changing the distribution @ . We also map 6 • 2 n ~ 3 words 
1-to-l, and finally we join two multiplicities #(vi),#(v2) by mapping vi,V2 onto the same w, the actual 
effect of this application of AS*. The new values @+ are thus @+(#(vi)) = @(#(vi)) - 1, @+(#(v 2 )) = 
@(#(v 2 )) - l,@ + (#(vi)+#(v 2 )) = @(#(vi)+#(v 2 )) + l, and @+(Jfc) = @(k) otherwise. In this way, 
we eventually arrive at a distribution @ as required by /. 

To achieve this, we permute values in between applications of AS*. In this (slow) way, we eventually 
get to the distribution of #(w) required by /. 

The final permutation %i maps the v with multiplicities #(v) > to the correct values w G Im(f). 

Since we always have two words v\ , V2 mapping to the same w under AS* , and also two words w\ , W2 
outside 7m(AS*), any n k G S2AA2" can be extended by one of the transpositions (vi,V2) or (w\,W2) to an 
equivalent n k G A2" . 

Concerning the "only if" part, already the first application of AS* would decrease the number of 
values below \Im (/) | . □ 

3.3 Examples 

The superscript W indicates the torus size. 

INC ^ 3 ) For n = 3, let w = v + 1 mod 8. This is an odd bijective function, and hence representable 
for this n = 3. 

MUL-BY-3 < 3 ) For n = 3, let w = 3 • v mod 8. Same as with INC. 

MUL-BY-2_( 3 ) For n = 3, let w = 2 ■ v mod 8. From 2-0 = 2-4 = mod 8, this function is not 
bijective, and hence not representable by ECA-57 for n = 3. 

INC( 4 ) For n = 4, let w = v+ 1 mod 16. As with n = 3, this is an odd bijective function. Contrary 
to the case n = 3, a representation by ECA-57 is not possible for n > 4. 

INC' ( 4 ) For n = 4, let w = v + 1 mod 15, 15 >->■ 15. This is an even bijective function, and thus 
representable. 

MUL-2-BY-2 W For n = 4, let v = a\b, < a,b < 3 and w = a ■ b. The range is given by the multiset 
{0 7 , 1,2 2 ,3 2 ,4,6 2 ,9}, where superscripts show the number of occurrences. The sum £ weF 4 L#(w)/2J = 
6 > 2"~ 3 = 2 is large enough (the range is sufficiently small thus) to allow shrinking by e.g. repeated 
application of the asynchronicity pattern "<=>>" and suitable permutations. Multiplication can thus 
be computed by ECA-57 through asynchronicity. How to do it exactly, is a more complicated case, see 
Open Problems. 

MUL-k-BY-k ^ 2 ^ : Zero appears 2-2^ — 1 times, and \ <a<b<2 — 1 yields a-b = b ■ a that is at 
least ( 2 *~ 1 ) pairs. Hence, we have £ wGF 4 [#(w)/2\ > 2 k - 1 + (2 k - 1) ■ (2 k - 2)/2 > 2 2 ^ 3 = 2"~ 3 (with 
k > 2). All these multiplications can therefore be computed by ECA-57, using asynchronicity. 
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Boolean and arithmetic functions on k bits, n = 2k: 

Let v = a\b with 0<a,b< 2k. Then / v ,i ( v ) = (0|a Vfc), /v,2(v) = {a Vb\a Vfe), / A ,i(v) = (0|aAfc), 
/a,2(v) = (qAfc|aAfc),/s,i(v) = (0|fle*Or/"e,2(v) = (a©%©&), A(v) = (o-fe mod 2 2 *), /_,(v) = 
(0| (a — b) mod 2 k ), can all be computed by ECA-57 under asynchronicity, for any k that is any even n: 

All these Boolean functions are commutative, a o b = b o a, and thus enough pairs (vi,V2) with 
/( Vl ) = /( V2 ) exist to have LL # (w)/2j > 2 2 *~ 3 . 

For < v < 2", let /neg(v) = (— v) (2's complement). This is an odd bijection with the two fixed 
points and 2 n_1 , and 2" _1 — 1 transpositions, hence computable by ECA-57 only for n = 3. 

/comp(v) = (v© 111 • • • 111) (l's complement), on the other hand, is an even bijection, computable 
for all n. 

Example 5 A detailed description of the calculation of INC^ 3 ' and MUL-BY-3 1 - 3 '. The left column 
indicates the temporal rule and partition of Z/3Z 

INC (3) 

000 001 010 011 100 101 110 111 

<><(1|0|2) on 1Q1 10 q n0 nl 01Q q 01 QQQ 

<»(i|2|o) 010 110 on 000 001 101 100 111 

<>=(i|o,2) 101 000 110 oil 100 010 111 001 

=><(o,i|2) oiO 111 100 110 011 001 000 101 

»<(o|i|2) ooi 010 011 100 101 110 111 000 
MUL-BY-3 (3) 

000 001 010 011 100 101 110 111 

«>(2|i|o) 10 i oiO 111 100 110 011 001 000 

><>(2|o|i) ooi 101 100 010 011 000 110 111 

=<>(2|o,i) 110 oil 010 111 000 101 001 100 

>=<(0|i,2) 101 10 o 001 010 110 000 111 011 

»<(o|i|2) ooo 011 110 001 100 111 010 101 

Table 5: INC (3) and MUL-BY-3 (3) in detail 



Further Research and Open Problems 

1. Give an algorithm to calculate the temporal sequence (ASi, AS2, . . . , AS^) for a function on F 2 directly 
from the function values, given e.g. as permutation on {0, 1, ... ,2 n — 1}, instead of searching through 
the full tree AS*. 

2. Consider temporal sequences that do not depend on the position, but on the rule to be applied, 
e.g. first update at all corresponding sites 000 >->■ p , then 110 h-> p 6 , then 010 >->■ p 2 etc. There are 8! = 
40320 such temporal rules, independent of n. 

How do we treat actions that already had their turn, but whose neighborhood only turns up later? 
Update immediately upon creation, never in this round, ...? 

This could mimic chemical reactions, e.g. in cell biology, DNA expression, where some reactions 
are faster than others, depending on their reaction rate constant k. 

3. As in 2., but associate a latency time with each temporal rule: As soon as the corresponding 
neighborhood pattern is created, wait for its latency time and then update according to the temporal rule. 
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4. What can we say about the alphabet {0, 1,2} instead of {0, 1}? 

There are now 3 3 « 2 43 ECA's to be considered. Since there are 3" — 2" +1 + 2 = 0(3") asynchronic- 
ities (Definition 1) and exactly 3" configurations on Z/nZ, an analogue of properties (o) to (iv) is now 
impossible due to lack of temporal rules. However, we may ask, how the number of configurations 
actually reached grows with n. Do we ever obtain the full diversity of 3" — 2" +1 + 2 results? 

Conclusion 

We have introduced temporal order via temporal rules as a means to diversify the behaviour of elementary 
cellular automata. 

In particular, ECA's with update rules 19, 41, and 57 are pattern universal for n < 12, achieving any 
desired pattern transduction vh^w, applying iteratedly a single temporal rule. We conjecture that they 
are indeed uniformly pattern-universal. 

ECA-57 produces any even (as permutation) bijective function on {0, 1}", for n < 10, and all non- 
bijective ones that join at least 2"~ 3 pairs of argument values. 

Temporal order is thus a third way to encode information and algorithms, after programs 
(ECA's) and data (initial configurations). 

This may have farreaching consequences, e.g. for modeling gene expression, since physico-biological 
processes seldomly achieve exact synchronicity. 
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Appendix A - ECA Families 

Each family (equivalence class under the symmetries 0/1 and L/R (c,_i o c,+i)) consists in up to 4 EC As with 
numbers EC A=abcdefgti2, 



ECAH^ aecgbfdh 2 , ECA &\ hgfedcba 2 , 


and ECA^i^^ 1 hdfbgcea 2 . 


(255) 


1 (127) 


2 (191 16 247) 


3 (63 17 119) 

-J \ Uy -1 / -1 -1 J 7 


4 (223), 


5 (95) 


6 (159 20 215) 


7 (31 21 87) 


8 (239 64 253) 


9 (1 1 1 65 125) 


10 (175 80 245) 

1 v l 1 / J \J\J T^y /* 


1 1 (47 81 117) 


12 (207 68 221) 


13 (79 69 93) 


14 C143 84 213) 


15 (85), 


18 (183) 


19 (55), 


22 (151), 


23, 


24 (23 1 66 1 89) 


25 (103 67 61) 


26 (167 82 181) 


27 (39 83 53) 


28 (199 70 157) 


29 (71) 


30 (135 86 149) 


32 (251), 


33 (123) 


34 H87 48 243) 


35 f 59 49 1 1 5) 


36 (219), 


37 (91) 


38 (155 52 211) 


40 (235 96 249), 


41 (107 97 121), 


42 (171 112 241) 

11/1 J. J. £^ I -11* 


43 (113), 


44 (203 100 217) 

1 1 \ £^\J *J -1 \J\J -1 / / j 


45 (75 101 89) 


46 (139 1 16 209) 


50 (179), 


51, 


54 (147), 


56 (227 98 185), 


57 (99), 


58 (163 114 177), 


60 (195 102 153), 


62(131 118 145), 


72 (237), 


73 (109), 


74 (173 88 229), 


76 (205 76 205), 


77, 


78 (141 92 197), 


90 (165 90 165), 


94(133), 


104 (233), 


105, 


106(169 120 225), 


108 (201), 


110(137 124 193), 


122(161), 


126 (129), 


128 (254), 


130(190 144246), 


132 (222), 


134(158 148 214), 


136 (238 192 252), 


138 (174 208 244), 


140 (206 196 220), 


142 (212), 


146(182), 


150, 


152 (230 194 188), 


154(166 210 180), 


156 (198), 


160 (250), 


162(186 176 242), 


164 (218), 


168 (234 224 248), 


170 (240), 


172 (202 228 216), 


178, 


184 (226), 


200 (236), 


204, 


232 



